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ABSTRACT 

Isolated black holes (IBHs) are not usually considered to be important astrophysical sources, 
since, even in the case of a high accretion rate, an accretion disk rarely can be formed due to 
the small angular momentum of the in-falling matter. Thus, such systems are not expected to 
^ c| feature thermal disk emission which makes the dominant contribution to the radiative output 

Q_f of binary systems harbouring a BH. Moreover, due to their relatively modest accretion rates, 

these objects are not conventionally treated as feasible jet sources. However, the large num- 
ber of IBHs in the Galaxy, estimated to be ~ 10 8 , implies a very high density of lCT 4 pc~ 3 
and an average distance between IBHs of ~ 10 pc. Our study shows that the magnetic flux, 
accumulated on the horizon of an IBH because of accretion of interstellar matter, allows the 
Blandford-Znajeck mechanism to be activated. Thus, electron-positron jets can be launched. 
We have performed 2D numerical modelling which allowed the jet power to be estimated. 
Their inferred properties make such jets a feasible electron accelerator which, in molecular 
clouds, allows electron energy to be boosted up to ~ 1 PeV. For the conditions expected in 
molecular clouds the radiative cooling time should be comparable to the escape time. Thus 
these sources can contribute both to the population of unidentified point-like sources and to the 
local cosmic ray (CR) electron spectrum. The impact of the generated electron CRs depends 
on the diffusion rate inside molecular clouds (MCs). If the diffusion regime in a MC is similar 
to Galactic diffusion, the produced electrons should rapidly escape the cloud and contribute to 
the Galactic CR population at very high energies > 100 TeV. However, due to the modest jet 
luminosity (at the level of ~ 10 35 erg s _1 ) and low filling factor of MC, these sources cannot 
make a significant contribution to the spectrum of cosmic ray electrons at lower energies. On 
the other hand, if the diffusion within MCs operates at a rate close to the Bohm limit, the CR 
electrons escaping from the source should be confined in the cloud, significantly contributing 
to the local density of CRs. The IC emission of these locally-generated CRs may explain the 
variety of gamma ray spectra detected from nearby MCs. 
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1 INTRODUCTION 

The Galactic population of isolated compact objects - rem- 
nan ts of massive star evolu tion - is very large: about 10 9 (see, 
e.g. ISartore& Treves! |2010| . and references therein). It should be 
dominated by neutron stars, but the number of isolated black 
holes (IBHs) is also non-ne gligible: their number is > 10 8 
jAgol & Kamionkowskil2002h . Taking into account this large num- 
ber, the local spatial density of IBH can be high, corresponding to 
a typical distance to the closest objects ~tens of parsec. Such es- 
timates are usually based on detailed calculations for isolated neu- 
tron stars and normalization based on the relative number of com- 



* E-mail: sergepolar@gmail.com 



pact objects. The local isolated neutron st ar spatial density is esti- 
mated to be ~ 3 x 1CT 5 pc~ 3 -3 x 1(T 4 pc 3 jPopov et al.l2005l : IOfekl 
l2009l : ISartore & TrevesboTol) . The relative number of BHs to neu- 
tron stars is about 1/3. Then, the local density of IBHs can be taken 
as ~ 1(T 5 - 1(T 4 pc" 3 . 

Historically, the first ideas about the detectability of IBH 
were related to a ccretion from the interstellar medium (ISM) 
dShvartsmanlll97ll) . Detailed studies concluded that ac creting IBH 
must emit mostly in X-rays and IR (see, for example, iFuiita et aU 
Il998l) . Calculations of the observability of the population of Galac- 
tic accreting IBHs were pe rformed several times (see, for example, 
IPodov & Prokhorovll 19981 and refe rences therein). One of the mos t 
detailed studies was prese nted by I Agol & KamionkowskH J2002h . 
Despite several attempts dChisholm et al. 2003 ). up to now no 



2 M.V. Barkov et al. 



firm candidat es for accreting IB Hs (or for accreting isolated neu- 
tron stars, see lTreves et al Hooch ar e known. Optimistica lly, future 
space missions may discover them dGrindlav et al .1120091) . 

The only possible detections of IBHs reported in recent years 
are related to microlensing. About 10 years ago s everal can- 
didates were reported by different researc h groups dAgol et al] 
120021; iMao et aljl2002l; iBennett et ai1l2002l) . Recent calculations 



1 Sartore & Treves 2010h predict that further compact objects will 
be detected via microlensing. 

To understand IBHs better, it is very important to seek ex- 
amples of these not related to microlensing events. There is a dis- 
covery potential for IBHs via mo delling of disrupted massive bi- 
naries dProkhorov & Popov! l2002h . or in a close examination of 
unidentified objects with peculiar properti es among different sur- 
veys. In this note we develop the idea of iBarkov & Khangulvanl 
d2012h who proposed that BHs with spheric al accretion — i.e. with- 
out fo rmation of an accretion disk (see also iKomissarov & Barkovl 
120091) — can be sources of relativistic jet s powered by rotating BH 
dRuffini&Wilsoij|l975l ; lLovelacdll976l) or. e.g.. th e Blandford- 
Znajek (BZ) mechanism iBlandford & Znai elj 1 19771) 1*1 These jets 
form conditions for efficient acceleration of non-thermal particles 
which can radiate energetic electrons, with quite low luminosity 
in X-ray and other energy bands. Strict limits on the energy re- 
lease rate in X-rays L x < 10~ 3 (in units of the Bondi-Hoyle ac- 
cretion rate with efficiency 0.1c 2 ) was achieved in the work of 
I 1 J 

IMotch & Pa kull (201 2D, which is complim entary to the scenario 
suggested bv lBarkov & Khangulvanl ( 120120 . 

In the following two sections we briefly describe the model. 
Then, in Sec|4]we present our results on the possible observational 
appearance of IBHs. Finally, we present some discussion and our 
conclusions. 



2 JETS FROM IBHS 

Recently, it was shown that the BZ mechanism can be activated in 
the case of d irect wind accretion on to a r otating BH in a close bi- 
nary system dBarkov & Khangulvanll2012T) . Thus, a relatively pow- 
erful jet, with kinetic luminosity of ~ 10 35 erg s -1 , can be launched 
in systems which do not feature an accretion disk. A similar sit- 
uation may arise in the case of accretion of the ISM on to an 
IBH. Accretion from the ISM has been studied f o r many decades 
startin g from the semi nal papers by lBondil ( [19521) : lBondi & Hovlel 
d 19441) (see a review in lEdgaj2004l). The case of IBHs was studied 



in dShvartsmanll97ll ; lBisnovatvi-Kogan & Ruzmaikinll974lll976h 



and later this appr oach was developed in 1 Beskin & Karpovj 2005 
iBeskin et alj2008h . 

In the simplest case, the accretion rate can be estimated as: 



M = XAn 



(GM) 2 p 



(1) 



Here M is the mass of an accretion source, p is the density of the 
surrounding medium (the ISM density in our case), and A is a di- 
mensionless coefficient < 1 . The velocity v must include the con- 
tribution from different types of motion. The most important is just 
a spatial velocity of the accreting object, u„. When this is small it is 
necessary to take into account the sound speed in the medium. Typ- 
ically, the sound speed in the ISM is ~ 10 km s _I . Spatial velocities 
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Punsfyl <1998bld)1999h : IPunslv et alj l2000h : iTorres et alj 
Maccarond i2005l) 



are larger, since progenitor stars have typical velocity distribution 
of Gaussian shape with cr ~ 30 km s , and the direction of the 
BH birth kick is not expected to be related to the direction of the 
progenitor star's spacial velocity. Thus, in this study we assume 
that the typical velocities of IBHs are about several tens of km per 
second (perhaps as high as hundred s km per second, sim ilar to the 
values typical for neutron stars, see lRepettoetai]|2012h . and con- 
sequently the additional contribution due to the ISM sound speed 
can be neglected. 

The rate of matter capture can be estimated from Eq.(QJ as 



M = 1 X 10"/t M, 1^-24 g s~ 



(2) 



here M\ = M/WM Q , u L5 = u/10 1 ' 5 km s _1 « u/31 km s" 1 , 
p_24 = p/10~ 24 g cm"" 3 are the mass of an IBH, the peculiar ve- 
locity of the IBH and the density of the ISM, respectively 0. Since 
IBHs are not expected to be powerful sources, we mostly focus 
on the most nearby lo cations dBocMtarev|[l99cl Isfeir et alj[l999l 
iMafz-Apelliiul 1200 lh . This allows the ISM density to be esti- 
mated. Namely, currently the solar system is located in the Lo- 
cal interstellar cloud, which is expected to have particle density 
of ~ 0.1 proton cm 4 , or equivalently p lic ~ 2 x 10~ 25 gcm~ 3 . 
The surrounding region, the so-called the Local Bubble, is char- 
acterized by a much smaller density, pib ~ 10~ 26 gcrrT 3 . Finally, 
at distances d mc ^ 150 pc from the solar system there is a num- 
ber of molecular clouds (MC) which may have rather high density 
Pmc ~ 10~ 21 gcirT 3 . Therefore, unless the BH proper velocity dif- 
fers significantly from v l 5 ~ 1, the expected accretion rates for 
the LIC, LB and MC are M hc 10" g s" 1 , M\i 



10 10 gs-\ and 



M mc a 10 15 g s~', respectively. 

In the cases that the accretion regime differs from Bondi- 
Hoyle accretion, the rate can be even smaller. In particular, for 
low velocity IBHs, u < 50 km s -1 , a temporary dis k can be 
formed due to interstellar turbulence dFuiitaetalJll998h . The ac- 
cretion rate in this case can be slightly smaller than the Bondi value 
dBlandford & Begelmanll 19991) due to the formation of a wind out- 
flow from the accretion disk. However, such disks are expected to 
be transient on a time-scale fbdt ~ Rb/voc, where Rb is the so-called 
Bondi radius: 



R B = 2GM/ui = 2.5 x 10 14 Miu7] cm 



and 



Ann 



Rb/Vc ~ 2GM/vi * lO 8 /^: 3 s. 



(3) 



(4) 



Thus, for realistic cases, this time-scale is very short, ~ 1 yr. There- 
fore in what follows we adopt the values obtained for the Bondi 

accretion regime. 

In the model bv lBarkov & Khangulvanl d2012t) a jet is driven 
by the BZ mechanism. The total power in the jet can be estimated 
by a simple numerical relationship 



L BZ ~ CMc 2 = 10 31 M H C-i erg s"\ 



(5) 



here C_i = C/0.1 is a non-dimensional parameter C_i ^ 1. The 
maximum value of C_i - 1 can be achieved if high enough mag- 
netic flux can be accumulated at the BH horizon. The required 
magnetic field flux can be ob tained through the following relation 
dBarkov & Komissarovll2008h : 

L BZ = 1.4 x VpfiayPvMferg s"\ (6) 

here = T/IO 17 is magnetic flux in units 10 17 G cm 2 . The 
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function/(fl) = cr(\ + Vl - a 2 } is a dimensionless function ac- 
counting for the BH rotatior0. Combining Eqs.l[5} and @ and 
adopting f(a) « 0.1 (or a ~ 0.5), one can derive the magnetic 
flux value on BH horizon, required for the optimal operation of the 
BZ mechanism: 



Tbz * 3 x 10 18 Af}f Afi Gem 2 



(7) 



The value of "P which an IBH can accumulate by Bondi accretion 

! 14^ism,-6<J ^rn . (8) 



is = 2R 2 3 B l&M or 



Therefore, in the case of accretion of the ISM matter on to an IBH 
the BZ mechanism can achieve its optimal performance provided 
by Eq. {5). However, we note that there are some uncertainties re- 
lated to role of the reconnection and magnetic flux escape from the 
horizon of an IBH. 

Thus, for the expected accretion rates, the IBH might be 
characterized by a jet with a luminosity of Lji ic ~ 10 31 ergs -1 , 
Ljjb ~ 10 30 ergs~' and Lj, mc ~ 10 35 ergs -1 for the cases of BHs 
located in LIC, LB and MC, respectively. Given the dependence of 
the jet luminosity on different factors, e.g., the accretion rate, on the 
magnetic field at the horizon, spin and velocity of the IBH, these 
luminosities can be, in fact, considered as upper limits. 



3 SPATIAL DENSITY OF IBHS 

Another important issue is related to the actual number of such 
IBHs in the ISM with different properties. For the above esti- 
mates there should be N lBli = 4 x (10 _1 - \Qr 2 )R\ located within 
R = 10R] pc from the solar system. For the expected size of the LIC 
of ~ lOpc, this yields a probability of 4% - 40% for a BH to be 
located there. With increasing R the number of BHs rises rapidly 
achieving 10 3 ~ 4 inside ~ 1 kpc. On the other hand, IBHs in high 
density regions are rare: they are at least two orders of magnitude 
less frequent than IBHs in the normal and low-density ISM. A typi- 
cal MC with M ~ 10 5 M Q has a volume ~ 10 4 pc 3 . I.e., for the IBH 
spatial density close to the Galactic plane ~ 10~ 4 pc~ 3 , we expect 
one IBH per cloud. On larger scales, the number of IBHs located 
in MCs might be very significant. Given, however, the presence of 
pulsars on this scale, the impact of IBHs is negligibly small. For 
example, the Vela pulsar is located at a distance of 290 pc, and is 
expected to be a powerful source of cosmic ray (CR) leptons with 
power co mparable to its spin -down luminosity, L 8 a = 7x 10 36 erg s _I 
(see e.g. iHinton et alj 1201 ll . and references therein). The closest 
IBH in a cloud is expected at a distance ~ 170 - 2 00 pc (see the list 
of near-by MCs, for example in lDame et al .1 19871) . Up to ~ 300 pc 
only ~ 5 IBHs are expected to be located in MCs. 



4 OBSERVATIONAL APPEARANCE OF IBHS JETS 
4.1 Particle acceleration 

Non-thermal processes in BH jets can be linked to shocks at the 
head of the jet. A system consisting of two shocks is expected to 
be formed: a forward shock propagating through the ISM and a 



3 Let us note that a BH can be significantly slowed down if it had accreted 
about 0.1 of its mass. In the case of IBHs using Eq.(2} it will last up to 
f sd ~ 0.1Mbh/M ~ 10 14 MiMi2 years, so we can neglect the IBH braking. 



relativistic reverse shock in the jet matter. These shocks are char- 
acterized by very different velocities. Properties of the non-thermal 
particles accelerated by these shocks may differ significantly there- 
fore. The jet power and medium density determines the key prop- 
erties of the forward shock. Namely, the jet ram pressure should be 
balanced by the inertia of the ISM. This yields the following rela- 
tion Lj = p mm dV/dt = p tam nR\c, where p mm ~ pv 2 and R$ are the 
jet ram pressure and the radius of the termination shock, respec- 
tively. The latter can be estimated as 



Rst 



3 x lO^i'^pIjfv^cm. 



(9) 



An upper limit on the forward shock magnetic field strength can be 
derived from pressure balance B^/Sk < p ram : 



15pi^wi.5j"G. 



(10) 



The reverse shock in the jet has a similar size and, in fact, its 
magnetic field strength should be comparable. Indeed, since jets 
launched by the BZ mechanism are magnetically dominated, the 
magnetic field strength can be estimated as 



/4Lj 1 
~c~~hll' 

This allows one to obtain the following relation: 

B,R S = W n Ml /2 v-;l /2 p l ^Cl?Gcm. 

For the shock size R$ determined by Eq.(|9), this yields a magnetic 
field of 



(11) 



(12) 



Sj = lOpi'^j/uyuG, 



(13) 



which is quite close to the strength expected at the forward shock 
Eq ,d lot- The main difference between these two shocks is related 
to the efficiency of the acceleration process. 

The acceleration time of non-thermal particles can be esti- 
mated as 



c 



(14) 



here R g is electron giro-radius, and 77 > 1 is a dimensionless pa- 
rameter, which in the case of diffusive shock acceleration has a 
value of j/jjff = 27r(c/u s h) 2 . In the case of the forward termina- 
tion shock, it is natural to assume that v s h m v. In the case of the 
reverse shock, which is expected be characterized by relativistic 
velocities, the value of ^-parameter is not defined in the framework 
of theoretical models. However, from the interpretation of the non- 
thermal emission detected from the Crab Nebula, one can infer that 
in the relativistic case rj < 100, thus in what follows we will assume 

>h* = 10. 

In the case of diffusive shock acceleration the maximum en- 
ergy of the non-thermal particles can be limited either by loss rate 
or by the confin ement requir ement. The latter is refereed to as the 
Hillas criterion dHillasll984T) . which for the Bohm diffusion regime 
(i.e., the diffusion coefficient is D ~ R s c) can be formulated as a 
limitation of the size of the acceleration site R., cc : 



R^>R g ^- 



(15) 



In the asymptotic case 77 = 1 so this condition results in an obvious 
requirement of R xc > R^. Since in the case of the jet termination 
Race ~ Rs, the electrons can be accelerated up to 



< SCifp^M^^GeV 



(16) 



This value gives ~ 2 GeV for the Local cloud, and ~ 0.5 GeV 
for a BH located in the Local Bubble. In the case of a BH in a 
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MC, the maximum energy achieved at the jet termination shock is 
~ 200 GeV. Therefore, the jet forward termination shock is not a 
plausible acceleration site for any reasonable conditions. 

Regarding the reverse shock, the Hillas Criterion gives the fol- 
lowing estimate: 

£max, < 30C 1 _[ 2 p%Ml / \l l2 T,- l/2 TeV . (17) 

This gives ~ 10 TeV for the Local cloud, and ~ 1 TeV for BH lo- 
cated in the Local Bubble. In the case of an IBH in a MC, the max- 
imum energy of particles accelerated at the jet termination shock is 
~ 1 PeV. 

4.2 Source emission 

The acceleration process is limited by the Hillas Criterion, there- 
fore the particles are not expected to loose all their energy in the 
source. Nevertheless, we compare the radiation cooling time to the 
particle escape time. The latter depends on the regime of diffusion 
inside the source. Obviously, this diffusion rate is different from the 
Galactic diffusion coefficient used to derive Eq. J27b and it cannot 
be determined from first principles. However, the X-ray properties 
of supernovae remnants apparently indicate that inside these shock 
related non- thermal sources, the diffusion regime is similar to the 
Bohm limit jAharonianir2004l) . We therefore assume this diffusion 
regime to be valid inside the source. Thus, the source escape time 
can be estimated as 

r bd = 5x 10 7 £ T >7.5Pi24 s ■ ( 18 > 

Moreover, we note that the assumption of Bohm diffusion inside 
the source is in fact closely related to the acceleration rate assumed 
above. 

The radiative cooling time is determined by the density of the 
corresponding target. In the case of the synchrotron and inverse 
Compton (IC) radiation mechanisms on the magnetic field e b = 
B 1 j%ix and photon £ ph energy densities: 

t = 3x 10 5 (£ ph /leVcm-V£ Ti ! v yr (19) 

and the synchrotron cooling 

f syn = 4 x 10 5 (S/5//G)- 2 £ Te ' v yr . (20) 

Equations [10] and Q~3] allow the synchrotron cooling time to be ob- 
tained, and the IC cooling time can be estimated based on the size 
of the termination shock. Jet kinetic energy is largely transformed 
to the thermal energy at the termination shock. The photon energy 
density inside the source can be estimated as £ ph ss L- s l2nR\c ~ 
3p_24Vj 5 eV cm 4 , which is comparable to the magnetic field en- 
ergy density: £ p h ~ s B ~ 5p_24f 2 5 eV cm 4 . Thus, assuming that 
IC scattering operates in the Thomson regime we can estimate the 
radiative cooling time as: 

f rad = 1.5 x lOVVrl^v s . (21) 

The lower energy part of the spectrum, E < 1 GeV, is dominated 
by bremsstrahlung losses 

f bs = 6x 10 7 p: 24 yr. (22) 

The non-thermal particles emit in the source a fraction of their 
energy, which depends on the ratio of the diffusion and radiative 
cooling times: 

^ (23) 

1 + 'bd/'rad 



where x is the fraction of the jet kinetic luminosity which is trans- 
ferred to the population of the non-thermal particles. The esti- 
mate Eq. d23l > depends strongly on the density of the surrounding 
medium: 

^ = 3xlO- 5 p!^7 2 5 . (24) 

>rad 

It can be seen that, for the case of the density of the medium ex- 
pected in MCs, i.e., p ~ 10 41 gem 4 , this estimate approaches 
unity. Thus, while in the cases of Local Bubble and Local inter- 
stellar cloud, jets emerging from IBHs are not expected to be plau- 
sible point sources, in the case of an IBH accreting in a MC, non- 
thermal particles can emit a significant part of their energy in the 
close vicinity of the source. One should expect similar fluxes emit- 
ted through the synchrotron and IC channels, since the energy den- 
sities of the target fields are similar. However, these components are 
located in different energy bands. The synchrotron emission should 
be emitted at an energy 

e syn = lOvuP^XeveV, (25) 

where £ TcV and p_ 2 i are the energy of the emitting particle in TeV 
units and density of the MC in 10 gem 4 units. For the assumed 
condition at the reverse shock electron energy can significantly ex- 
ceed 1 TeV (see Eq. [T6l >, thus the synchrotron emission can cover 
quite a broad range from radio to soft gamma-rays. 

The main target for the IC scattering is provided by the 
shocked MC material, which is expected to be heated to ~ 50 K. 
In the Thomson limit the energy of the up-scattered photon can be 
estimated as: 

e ic = 20£ 2 W GeV . (26) 

The transition to the Klein-Nishina regime occurs for electrons 
of energy > 20 TeV, i.e. in the band which is slightly above the 
range of optimal sensitivity of modern Cherenkov detectors. Thus, 
in what follows we do not discuss this energy range. It is also wor- 
thy to note that the spectral shape is expected to be flat with photon 
index ~ 2, without any significant breaks. Indeed, since the domi- 
nant cooling processes - synchrotron and IC - and escape time have 
the same energy dependence, the spectral breaks can be induced ei- 
ther by the transition to bremsstrahlung cooling dominance, or by 
the age of the source. The later is not relevant for sources with age 
> 10 4 yr, and the bremsstrahlung break should appear in the syn- 
chrotron spectrum at energy of < 1 GHz, i.e. in the range, which 
can hardly be probed with observations. 

However, if the diffusion deviates significantly from the Bohm 
regime, this may lead to modification of the spectral shape, and a 
break in the electron spectrum should appear at energy correspond- 
ing to fesc = Also a weak Klein-Nishina harden ing (see e.g. 
iKhangulvan & Aharoniai"j|2005l : iModerski et alj|2005h is expected 
to appear in the X-ray energy band e syn ~ 10 keV. 



4.3 Particle cooling and diffusion in ISM 

Since even in the case of IBHs located in MCs a significant fraction 
of the non-thermal power is able to escape from the source, below 
we study whether these sources can produce a detectable contribu- 
tion to the diffuse background. However, the total power of these 
sources per kpc 3 is ~ 10 36 ergs 4 kpc 4 ; this is approximately 2 or- 
ders of magnitude below the rate req uired to power the CR electron 
spectrum (see e.g. lAharonianl 120040 . Therefore, these sources can 
give an important contribution only in the high energy part of the 
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spectrum, where more distant sources are heavily suppressed due 
to particle energy losses. 

Following previou s studies l lPorter et al.l 120081 : IStrong et al.l 
120001 : iDwek et al]| 19971) we can say that the Galactic background 
radiation field in the solar vicinity is determined by 3 main 
components: 1) cosmic microwave background (CMB) (u cmb = 
0.25 eVcirT 3 ); 2) far infra-red (FIR) radiation of a dust (w flr = 
0.3eVcirT 3 ) and 3) visible light radiation from stars (u op = 
0.4eVcirr 3 ). We should point out that for optical radiation, lepton 
with E > ITeV interact in the Klein-Nishina regime, this signif- 
icantly suppressing the interaction cross-section. The energy den- 
sity of the magnetic field ne ar the Sun is u B = B 2 /Sn = 0.6 eV cm 3 
(see, e.g., lBochkarevlll990l) . The total energy density can be deter- 
mined as m 1+ b = "fir + "cmb + "b + «op ~ 1.2—1.5 eVcm -3 . 

The cooling time scales (Eqs. d 191221 should be compared to 
the diffusion time 

fair = W 5 R 2 2 5 E T %yr. (27) 

Here R2.5 = R/310 pc and 6 ~ 0.5. This diffusion time is normal- 
ized to the distance to the closest young energetic pulsar (the Vela 
pulsar) which is ex pected to be a powerful source of CR electrons 
dHinton et al 1 l201 lh . A comparison between f; c and ?dif allows one 
to determine the energy range for which the nearby sources can 
give a more important contribution than pulsars. Namely, a cooling 
cut-off is expected at energy 

E cat a 5R^f ^ 6) TeV. (28) 

In the case of IBHs (R=170 pc, see Eq.[3} R 2 ,s = 0.56 and we get 
an exponential cut-off near 50 TeV. Taking into account that the 
relation between pulsar and IBH power is £ ss 100, we can estimate 
the energy range where IBHs should dominate the spectra of CR 
electrons and positrons as £"b, ~ E cut log(^) 1/(1_l5) as 100 TeV. On 
the other hand, these energies are close to the limiting value for 
IBHs. If however multi- TeV CR electrons are d etected (e.g., with 
Cherenkov detectors, see lAharonian et al]|2008l) . the IBH scenario 
could provide a feasible explanation. Future experiments like CTA 
can shed light on this problem. 

Regarding a possible contribution to CR positrons, IBHs can 
supply particles in the high energy band (GeV and TeV). However, 
their contribution is expected to be significantly exceeded by the 
contribution from the nearby pulsars, e.g. from the Vela pulsar. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we examine IBHs as possible jet sources. We find 
that the non-thermal efficiency of these sources depends strongly 
on the density of the surrounding medium. IBH jets launched in 
high density environments are characterized by a larger power and 
a stronger magnetic field. The latter feature appears to be very 
important since (1) it allows non-thermal particles to be acceler- 
ated up to VHE; and (2) the ratio of the escape time to radiative 
cooling time decreases with increase of the magnetic field. The 
key source parameters (jet luminosity, maximum energy of non- 
thermal particles and non-thermal luminosity) scale with density 
as: Lj oc p, £ max oc p 1 ' 2 and L p h oc p 5/2 . The strong dependence of 
the non- thermal power, L ph , on the medium density is related to the 
change of the ratio of radiative to escape losses and, in the case of 
high enough density, converges to the fast cooling regime, where 
L ph OC Lj oc p. 

More specifically, if the density of the surrounding medium is 
high, p > 10~ 21 gem -3 , the radiation cooling time is shorter than 



the escape time (if the Bohm diffusion regime is applies within the 
source). Thus, the jet reverse shock can turn into an efficient emit- 
ter of non-thermal radiation. On the other hand, if the density is 
lower, p < 10~ 21 gcnT 3 , the IBH should act as a source of elec- 
tron CR, interestingly, for conventional values of the MC density 
(p ~ 10~ 21 gem -3 ), the non-thermal power of these sources should 
be shared in comparable fractions between the radiative cooling 
inside the source and escape. Finally, for the IBHs accreting in 
low density environments p <k 10~ 21 gem -3 , the accretion rate and 
maximum available energy of non-thermal particles appears to be 
too low to make any detectable contribution to the CR spectrum. 
Moreover, as fast particles escape these sources, then they should 
not be considered as feasible point-like emitters. 

As shown above, in the case of accretion in high density 
medium, the jets launched by IBHs, may emit a significant frac- 
tion of the accreted energy through synchrotron (from radio to soft 
gamma-rays) and IC channels (up to very high energies) within 
10 16 cm from the jet termination point. The luminosities radiated 
via these two mechanisms should be similar, since the energy den- 
sities of the target fields are expected to be comparable. The dom- 
inant target photon field for the IC is expected to be provided by 
shocked MC matter, which should be heated to T ~ 50 K. Because 
of this relatively low temperature of the target field, the Klein- 
Nishina effect should not affect the shape of the IC spectrum in 
the range relevant for the modern gamma-ray instruments. How- 
ever, the Klein-Nishina effect may lead to a weak hardening of the 
synchrotron X-ray compone nt (see, e.g jKhangulvan & Aharonianl 
I2OO5I : iModerski et al.ll2005l) . Thus, the spectra are expected to be 
featureless, with a photon index of ~ 2 over a broad energy range. 
Any stronger spectral features in the spectrum would be attributed 
to a deviation of the escape rate from the Bohm diffusion regime. 
It is also important to note that one can expect a certain variability 
of the non-ther mal flux caused by the inho mogeneous structure of 
MCs (see, e.g.. iKaplan & Pikel'Nej[l979h . For a typical velocity 
of 30kms -1 , these inhomogeneities should lead to variability on a 
year time scale. 

For the expected non-thermal luminosity of < 10 35 ergs~', 
such sources should be detectable with modern X- and gamma- 
ray instruments up to a distance of ~ 1 kpc. Within this distance 
~ 10 IBHs should be located in MCs. This may provide a feasible 
interpretation for some unidentified gamma-ray sources. 

Regarding the impact of the non-thermal particles escaped 
from the shock region generated by IBHs in MCs, there are cer- 
tain uncertainties related to the diffusion of these particles. Namely, 
the propagation of such particles may proceed very differently 
within and outside the MC. If the diffusion coefficient in the MC 
is comparable to the Galactic one, then the high energy parti- 
cles should rapidly escape from the MC (the MC crossing time is 
~ 300 yr). However, several uncertainties do not allow us to claim 
IBHs to be the best candidates for production of very high cos- 
mic ray electrons. The distribution of spins of IBHs is unknown 
and if a < 0.3 jets will be weak or simply not launch ed at all 
dKomissarov & Ba rkov 2009; Ba rkov & Komissarovl201oT) . Never- 
theless, we have checked the contribution from the nearby pulsars 
up to the distance 300 pc using ATNF (The Australia tele scope na- 
tional facility) pulsar catalogue (Manchester et al. 2005). We can 
conclude that the contribution in the CR electron spectrum from 
close pulsars is about a factor 10 smaller than expected contribution 
from IBHs. Although we can miss some pulsars if their beamed ra- 
diation does not cross the Earth, such pulsars can contribute to the 
population of electrons near the Sun. So far it is unlikely that IBHs 
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are the main source of electrons in CR, except hypothetical parti- 
cles with E > 100 TeV. 

However, the physical scenario should be quite different if 
in MCs the diffusion operates slower. For example, for the Bohm 
diffusion regime, the particles confinement in MC should be very 
long, ~ 10 7 yr. Thus, the non-thermal particles injected by a IBH 
should produce a diffusive emission component associated with the 
host MC. The luminosity of this component should be compara- 
ble to the non-thermal power of the source, i.e. < 10 35 erg s , so it 
might give a contribution compar able to the fluxes produced by CR 
intera cting with the MC material jAha ronian 200ll; ICasanova et"al] 
|2010|) . Importantly, the spectra produce by CR protons interacting 
with the MC should be quite different from the IC spectra gener- 
ated by the local population of CR electrons. Thus, this scenario 
can contribute to our understanding of the variety of the GeV spec- 
tra obtained from different MCs with Fenni/LAT (Yang and Aha- 
ronian, in preparation). 
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